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Abstract 

We extend the reduction group method to the Lax-Darboux schemes associated with nonlinear 
Schrodinger type equations. We consider all possible finite reduction groups and construct correspond¬ 
ing Lax operators, Darboux transformations, hierarchies of integrable differential-difference equations, 
integrable partial difference systems and associated scalar partial difference equations. 


1 Introduction 

In the theory of integrable systems the connections between partial differential equations, differential- 
difference and partial difference systems are well known. A clear and seminal account of these connections 
can be found in 0, m- They can be formulated in the frame of a Lax-Darboux scheme, where 

• the Lax structure (Lax representation, also known as zero curvature representation) is associated 
with partial differential equations (PDFs) and their symmetries [261 IT]: 

• Darboux transformations, which are automorphisms of the Lax structure, lead to Backlund trans¬ 
formations which can be regarded as integrable differential-difference equations (DAEs) |17l 1161 17]: 

• Bianchi permutability of the Darboux transformations yields integrable partial difference equations 
(PAEs) whose symmetries are the former DAEs [281 1^. 

In this paper we extend the reduction group method m to Lax-Darboux schemes for nonlinear 
Schrodinger type equations. More precisely, we study Lax operators of the form 

C = + U{p,q;X), (1) 

where the 2x2 matrix U belongs to the Lie algebra s[ 2 (C(A)). Matrix U{p,q]X) depends implicitly 
on X through two potentials p, q, and is a rational function in the spectral parameter A. Imposing 
the invariance of operator £ under the action of a reduction group, which is a finite subgroup of the 
group of automorphisms of s[ 2 (C(A)), we construct systematically the Lax operators corresponding to 
deep reductions. In this case there is a complete classification of finite reduction groups PEano] and 
corresponding reduced Lax operators [9]. Namely, in the sl 2 (C(A)) case there are only five distinct cases: 


(i) the trivial reduction group (no reductions); 

(ii) Z 2 group with a degenerate orbit; 
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(iii) Z 2 group with a generic orbit; 

(iv) Z2 X Z2 group with a degenerate orbit; 

(v) Z 2 X Z 2 group with a generic orbit. 

In the cases (i)-(iv) we construct an invariant Lax operator, a corresponding PDE, invariant Darboux 
transformations, corresponding integrable DAEs and PAEs. The simplest case (i) has been studied in 
detail in [7j. We present it here for completeness, in order to illustrate all elements of the corresponding 
Lax-Darboux scheme, such as dressing chains (also known as Backlund transformations) and their first 
integrals; to give a detailed derivation of associated integrable PAEs, and to discuss possible initial-value 
problems for these PAEs. The case (v) can be studied by the methods presented in the paper but leads 
to cumbersome expressions, and we have decided to omit it in order to keep our results presentable. 

Darboux transformations are automorphisms of the Lax structure and discrete symmetries of the cor¬ 
responding PDEs. With each Darboux transformation we associate an infinite lattice and a map. If there 
are two Darboux transformations, then the condition of their commutativity (the Bianchi permutability) 
yields an integrable system of PAEs. 

Although the theory of Darboux transformations is rather well developed and has a long history, 
there are a few important problems which require further research. One of the problems is to give 
a complete description of all possible Darboux transformations for a given Lax operator. In the case 
of the Schrodinger operator the solution is known; there is one Darboux transformation (depending 
on a parameter) and any other Darboux transformation can be represented as a composition of such 
transformations and their inverses for a certain choice of the parameters [5]. However the description of 
all possible Darboux transformations associated with a given Lax operator is still an open problem. 

The paper is organised as follows. In the following section, we introduce our notation and give the 
general scheme of these considerations. In the next four sections we consider the Lax operators related to 
the nonlinear Schrodinger equation (Section [3|), and operators derived from the reduction group method, 
pniE], namely Z 2 reduction (Sections H] and [5]) and dihedral group reduction (Section [6]). 

2 Lax-Darboux scheme 

In this section, we explain our terminology by describing the Lax-Darboux scheme. We present the 
class of Lax operators under consideration and discuss our general assumptions for the construction of 
Darboux matrices. Moreover, we introduce the notation we use throughout the paper. 

With the single term Lax-Darboux scheme we describe several structures which are related to each 
other and all of them are related to integrability. To be more precise, the Lax-Darboux scheme incor¬ 
porates Lax operators, corresponding Darboux matrices and Darboux transformations, as well as the 
Bianchi permutability of the latter transformations. 

• Lax operators are linear operators of the form L = -|- U, where the N x N matrix U is an 

element of a specific Lie algebra. As it was described in the previous section, in this paper we 
consider only the case where U{p, g; A) is a 2 x 2 matrix belonging to the Lie algebra s[ 2 (C(A)), and 
its dependence on the continuous variable x is implicit through the potentials p and q. 

• Darboux transformations S are automorphisms of the Lax operator C. They map T to T by 
updating potentials p and q. In other words, 

S'.C^C, where C = DxU{p,q-, X), C = DxU{p,q-, X), 
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with p, q denoting the updated potentials. 

Darboux transformations consist of Darboux matrices M along with corresponding dressing chains 
or Backlund transformations. 

• A Darboux matrix M maps a fundamental solution of the equation /l('I') = 0 to a fundamental 
solution 'h of /l('l') = 0 according to 'k = M^. In general, matrix M is invertible and depends on 
p, q, their updates p, q, the spectral parameter A, and some auxiliary functions. 

• Dressing chains are sets of differential equations relating the potentials and the auxiliary functions 
involved in L and L. They can be regarded as integrable systems of DAEs. This follows from 
the interpretation of the corresponding Darboux transformation as defining a shift on the lattice 
according to the sequence 



• If the Lax operator admits two commuting Darboux transformations S and T, then they define a 
two-dimensional lattice for which we adopt the multi-index notation {pij,qij) = {p,q), where 

i,j E Z. This interpretation allows us to derive systems of integrable PAEs by considering the 
Bianchi permutability of the corresponding transformations. 

In order to implement the above scheme, firstly we construct Darboux transformation S. From the 
definition of Darboux matrix M follows that 

MCM-^ = Z, (2) 

or, denoting the updated potentials with pio, qia and matrix U{pio,qio‘,X) with f7io, we can rewrite 
equation ([2]) explicitly as 

D^M + UioM - MU = 0. (3) 

For a given Lax operator C, the above equation can be used to determine M, as well as the corresponding 
dressing chain. Moreover, since matrices U and Dio are traceless, it follows from Abel’s theorem that 
the determinant of M is a first integral of the dressing chain. For the Lax operators we consider here, 
it is natural to assume that matrix M depends rationally on the spectral parameter A, and inherits the 
reduction group symmetries of the corresponding operator C. 

The interpretation of the Darboux transformation S as defining a lattice direction allows us to think 
the updated potentials in £ as shifts of the original ones p, q in that particular lattice direction. In this 
semi-discrete setting, the corresponding dressing chain can be seen as an integrable differential-difference 
equation da IS] deriving from the compatibility condition of the Lax-Darboux pair (also referred to as 
semi-discrete Lax pair) 


= -U{p,q;X)'i>, Tio = M(p, g,pio, gio; A)T . 

In this discrete interpretation, the Bianchi permutability of two different Darboux transformations 
yields an integrable system of PAEs in two discrete variables. Employing the standard notation for 
difference equations, we denote the two discrete variables with n and m, and interpret S and T as the 
corresponding shift operators defined by 

{h{n, m)) = h{n + i,m + j) = hij . 
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In particular, when i = j = 0, we will omit the index “00”, i.e. h = h{n,m). 
Now the shift operators <S and T act on a fundamental solution 'h as 


5 : 1-^ ^'lo = M{p,q,pio,qio-,f]X)^ = M 

T : ^01 = K{p,q,poi,qoi-,g-X)^ = K'^, 

where M and K are the corresponding Darboux matrices with / and g denoting any auxiliary (vector) 
functions. The Bianchi permutability of ([1]) according to Figure [T] allows us to compute Tn in two 
different ways. This yields the consistency condition 

T{M) K - S{K) M = (), (5) 

which is nothing else but the compatibility condition of the Darboux pair (also referred to as fully discrete 
Lax pair) 

5(T) = M{p,q,pio,qio;f-,X)^, T(T) = K{p,q,poi,qor, g] X)^ . 


^01 5 ^11 



Figure 1: Bianchi commuting diagram 

The resulting condition ([5]) yields a set of polynomial equations for p, q, /, g and their shifts. This 
set may have two branches of solutions. One of them leads to a trivial system, cf. (HID below, whereas 
the other branch yields a non-trivial integrable system of partial difference equations. Symmetries and 
first integrals for the non-trivial system follow from the dressing chain and the first integrals of the 
corresponding Darboux transformations. 

For some of these discrete systems, we employ first integrals and conservation laws to reduce the 
number of dependent variables and derive integrable scalar equations of Toda type. The form of these 
systems allows us to formulate a Cauchy problem on a single or a double staircase. 

In our derivations, we hnd more than one Darboux transformation for each Lax operator we consider. 
We would like to emphasise here that the interpretation of any pair of Darboux matrices as a discrete Lax 
pair as described above does not always lead to a non-trivial discrete system. In the following sections 
we present only the pairs of Darboux matrices which lead to genuinely non-trivial discrete integrable 
systems. 


3 Nonlinear Schrodinger equation 

In order to illustrate our approach, we consider a well known operator 

£ = Dj; + U{p,q-,X) = D^ + Xa 3 + 0-3 = diag(l,- 1 ), ( 6 ) 
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which is the spatial part of the Lax pair for the nonlinear Schrodinger equation [32] 


Pt = Pxx + 4:p^ q, qt =-qxx - ^Pq^■ (7) 

It is straightforward to verify that the constant matrix 

M=(“ a/3^0, 

is a Darboux transformation for this operator corresponding to the scaling symmetry of ([7]). 

pio = a,d~^p, gio = I3a~^q. 


The simplest A-dependent Darboux matrix one may consider is 


M = AMi + Mo. (8) 

Substituting ([ 6 |) and ([ 8 |) into the compatibility condition (l3|), the coefficient of A^ implies that matrix 
Ml must be diagonal. Additionally, from the diagonal part of the coefficient of A we conclude that Mi 
must be constant. Hence, Mi = diag(ci,C 2 ). We could choose either ci = 1, C 2 = 0, or ci = 0, C 2 = 1 
or Cl = C 2 = 1. Since the first two choices are gauge equivalent and the third one can be given as a 
composition of two suitable Darboux matrices with ciC 2 = 0, we choose 


Ml = 


1 0 
0 0 


Moreover, the off-diagonal part of the coefficient of A implies that the (2, 2) element of Mq is constant 
and, hence, we have to consider two distinct cases. 

The first case corresponds to 


M = A 



f a\ 

6 0 y ’ 


i.e. the (2,2) entry of Mq is zero. In this case, equation ([3]) is equivalent to the system 


a = p, b = qio, fx = 2{aq - bpio), ax = 2fp, b^ =-2 f qiQ. 


(9) 


The first two equations determine functions a and b, while the last two implies that pqio 
a non-zero constant (since det M = 7 0). Without any loss of generality we can set 7 = 

have 


1 

Q'lo — - j 
P 


pio=p (pq- ^fx 


p^ 

2p 


- 7 , where 7 is 
1 and thus we 

( 10 ) 


Finally, the Darboux matrix is given by 


M(PJ) = X[1 ") + 



(11) 


and the dressing chain (the Backlund transformation (llOp i can be rewritten in the form of the Toda 
lattice in a new variable (j) = logp 

(j)^^ = 4e'>-Wi,o _ 4e<Aio-<A. 
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In this case the Darboux transformation {p,q) —>■ {pio,qio) is explicit 


pio=p[ pq 


^ (Px 


4 \ p 


910 = 


P 


Alternatively, we can choose the (2, 2) element of Mq to be non zero and, without loss of generality 
set it to 1, i.e. 

/ a 

. b 1 

Now, it follows from ([3]) that 


Mo = 


a = p, b = qio, 

dxf = Hpq-pioqio), (12a) 

dxP = 2{pf - pio), dxqio = 2{q - qio f). (12b) 


A first integral of the above system is provided by the determinant of M, det M = A + / — pqiQ 

dxif -pqio) = 0- (13) 

Hence, matrix M has the following form 

= + (14) 

and (|12p is the corresponding dressing chain. 


3.1 Derivation of discrete systems 

Having derived two Darboux matrices for operator ([6]), we focus on the generic one given in (|14l) and 
consider the following Darboux pair 

I'm = M(p,gio,/)4', I'D! = M(p,goi,5')^, 


which explicitly reads as follows. 

The compatibility condition of (fT^ results to 


4'ni = A 


1 0 
0 0 


+ 


9 P 
901 1 


T. 


foi- f - (910 - 9 ) = 0 , 

/oi 9 - / 910 - Pio 910 + Poi 901 = 0 , 
P (/oi - 910 ) - Pw +P 01 = 0 , 

911 (/ - 9 ) - 901 + 910 = 0 . 


(15) 


(16a) 

(16b) 

(16c) 

(16d) 


This system can be solved either for (poi> 9oi>/oi j 9 ) or for (pio, 9 i 0 )/> 9io)- It has two branches of 
solutions. A trivial one 

Pw = Poi, 910 = 901) f = 9, 910 = /oi) (17) 
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corresponds to M{p,qio, f) = M{p,qQi,g), and a non-trivial solution given by 


_ qiop^ + (510 - f)p + pio 

Pro 5ii^ + (/-510)911+ 510 

(18a) 

POl 1 

l+pqii 

9oi .. 1 

l+pqii 

f qii (pio + pgio) + f - pqio 

Joi — , , 

l+pqii 

911 {pf - Pio) + 510 + pqio 

: 9 1 1 

l+pqn 

(18b) 


Some properties of the above system follow immediately from the derivation of the corresponding Darboux 
transformations. First of all, it admits two first integrals, cf. relation (jl3[) . namely 

{T - 1) if - pqio) = 0 and (<S - 1 ) (5 - pgoi) = 0 . (19) 

We can interpret functions / and g as being given on the edges of the quadrilateral where system (jl 8 jl is 
defined, and, consequently, consider system (fT 8 |) as a vertex-bond system m- System (|18l) admits the 
conservation law 

{T-l)f = {S-l)g 

which is the first equation in ()16ll . 

Moreover, relations (m imply that system (jl 8 p admits one generalised symmetry generated by the 
differential-difference equations 

dxP = fp- pio = gp- PQi, 

dxq = q -10 - f-m = qo-i - 50-W, ( 20 ) 

dxf = pq- pioqio, d^g = pq - pQiqQi. 


Evolution is determined by system il8t 


Evolution is determined by system i2H 




Figure 2: Initial value problem and direction of evolution 

Our choice to solve system (fT 6 ]l for poi) 9 oi) /oi and g is motivated by the initial value problem related 
to system (fT 8 ]l . Suppose that initial values for p and q are given at the vertices along the solid staircase 
as shown in Figure [2j Functions / and g are given on the edges of this initial value conhguration in a 
consistent way with the first integrals (1221) . In particular, horizontal edges carry the initial values of / 
and vertical edges the corresponding ones of g. With these initial conditions, the values of p and q can 
be uniquely determined at every vertex of the lattice, while / and g on the corresponding edges. This is 
obvious from the rational expressions (jlSp dehning the evolution above the staircase, cf. Figure [2j For 
the evolution below the staircase, one has to use 


_ qoip^ + if 01 - g)p + poi 

_ Poi 911 ^ + (5 - /oi) 911 + 901 

(21a) 

PlO 1 . 5 

1 +P911 

Q^io I ? 

1 +P911 

911 {Poi + pf 01) + 5 - pqoi 

510 - , 

1 +P911 

r _ 911 {pg - Poi) + foi + 5901 

1 J I ? 

l+pqil 

( 21 b) 
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which uniquely defines the evolution below the staircase as indicated in Figure [5J 

We could consider more general initial value configurations of staircases of lengths and (.2 in the 
n and m lattice direction, respectively. Such initial value problems are consistent with evolutions (jl8p . 
(|2T|) determining the values of all fields uniquely at every vertex and edge of the lattice. 

It follows from (fl^ that 


f-pqio = a(n) and g-pqoi = /3(m) 


( 22 ) 


and we can use these relations to eliminate / and g from (j21h . This results to a non-autonomous partial 
difference system for p and q only 

a{n) - I3{m) a{n) - I3{m) 

Poi = pio - ^r—zrz —9oi = qio H- ^ — qn- (23) 


1+pqii 


1+pqii 


Symmetries of this system can be derived directly from corresponding symmetries of system (|18l) by 
taking into account (1^ . In particular, it follows from (1^ that 

dxP = 2{p^ qio + a{n)p - pio), dxq = 2{q-io - p-io q^ - a{n - l)q) 


is a symmetry of 


3.1.1 Derivation of the discrete Toda equation 

Returning now to the construction of a discrete Lax pair, we employ matrix M{p, /), given in (jlip . and 
matrix M{p,qQi,g), in (fTT|) . That is, we consider the following system 


'I'lo — ( A 


1 0 
0 0 


+ 1 


/ P 


^ * 0 0 ' + 


g p 
qoi 1 




The compatibility condition of the above system implies that 

1 . ^ , 901 

P =—, 5 = a(m)H-, 

910 910 


as well as 


r 9oi 910 X r 911 910 , / N 

/ =- \-a{m), /oi =- \-a{m). 

910 911 920 911 


From the consistency of the latter equations and setting q = exp{—w-i-i), we derive the fully discrete 
Toda equation [miia] 




(24) 


along with its generalised symmetry 


dxw = - aim). 


Moreover, a conserved form of Toda equation is 

(r - 1) (e^o.-i-w-io _ ^ = [S- i)e^o,-i-w-io_ 

It is worth noting that a staircase initial value problem for the Toda equation (1241) involves the points 
Wi^-i and Wi-i-i, i.e. a staircase which is the reflection of the one shown in Figure [2] with respect to a 
vertical or horizontal line of the discrete plane. 
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4 Z 2 reduction group: Degenerate orbit 

Let us now consider an operator £(A) which is invariant under the transformation 


'Si(A) : C{X) —>■ asC{—X)a^. (25) 

The above involution generates the reduction group m which is isomorphic to the Z 2 group. The 
invariant operator corresponding to this orbit can be taken in the form 

C = D, + X^a3 + X(^^^ (26) 

and it is the spatial part of the Lax pair for the derivative nonlinear Schrodinger equation m 

Pt=Pxx + 4:{p^q)x, qt =-Qxx + 4:{pq‘^)x- (27) 

It can be easily verified that the constant matrix 

(28) 

is a Darboux matrix for operator (12611 corresponding to the scaling symmetry of system (12711 . 

pio = af3~^p, qio = Pa~^q. 


Considering Darboux matrix M with the same symmetry, i.e. M(A) = a^M{—X)a 2 „ we find after 
some analysis that the simplest A-dependent Darboux matrix can be written in the form 

M = A^Ma + XMi + Mq, 


where matrices M 2 and Mq are diagonal and matrix Mi is off-diagonal. Additionally, from the compati¬ 
bility condition ([3]) follows that Mq is a constant matrix. Moreover, following an argument similar to the 
one we used in the previous section, we consider only the case rank(M 2 ) = 1. Hence, summarizing the 
above analysis, we choose 

o)- "■ = (" 0) "»=(o "J- 


With these choices, equation ([3]1 firstly determines functions a, b in terms of /, p and qiq. In particular 
we find that 

a = fp, b = fqiQ. 

In terms of these relations, Darboux matrix becomes 


M{p,qio,f]Ci,C2) = A^ 


/ 0 
0 0 


+ A 


0 fp 

fqio 0 




Cl 0 A 

0 C 2 J ’ 


(29) 


and we derive the Backlund transformation 

dxP = 2p {pio qiQ-pq) -2 ^ ^ 

dxf = ‘2f {pq-pioqio) ■ 


2qio {pioqio -pq) -2 


ciqio -C 2 q 

f 


(30) 
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A first integral of the above system, which also guarantees that the determinant of matrix (I29p is inde¬ 
pendent of x, is given by 

dx {fpqio - C2f) = 0. (31) 

It is apparent that if constants ci, C 2 are not zero, we can always set them to 1 by composing Darboux 
matrix (I29p with an appropriate Darboux matrix ()28[) . Hence, we can impose without loss of generality 
that these constants are either 0 or 1. There are two particular sets of values for these constants at which 
differential-difference equations (I30p can be brought to polynomial form. 

1. First we consider the case when ci = C 2 = 0. It follows from equations (|30l) that f = 1/p and 
Qio = P) in view of which matrix M degenerates to 

M(p) = v(‘/f’ “)+a(; J). (32) 

The corresponding Backlund transformation becomes 

qio = P, dxP = 2 (pio - q) (33) 

and the first integral (I3ip holds identically. The resulting differential-difference equations (|33p are 
the modified Volterra chain. 


2. When ci = 1 and C 2 = 0, the Darboux matrix becomes 


Mip,qioJ) = 


/ 0 
0 0 


+ A 


0 fp 

fqio 0 




1 0 
0 0 


(34) 


the Backlund transformation simplifies to 


dxP = 2p {pioqio - pq) + y, d^qio = 2qio (pio^io - pq) - d^f = 2f{pq - piogio), (35) 

and the first integral (|3ip becomes 

dx{fpqio)=0- (36) 

In the context of differential-difference equations, if we make the point transformation 

p = u^, q = 


and subsequently, using the first integral (l36P . set 

j u V = I 

system ()35p can be written in a polynomial form as 

a /■ 2 2 2 2 \ \ 2 o /22 22\-.-2 

OxU = U[UiqV — U U_io) ±U V, OxV = V{UiqV — U U_;^g) =F UV . 


/ = —, 
uv 
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4.1 Derivation of discrete systems 

Now we consider the difference Lax pair 


= M(p,gio,/;ci,C2)4', ^'oi = , 


(37) 


where matrix M is given in (|29p and at least one of the constants ci, C 2 is different from 0. It follows 
from the above system that 


/5io - ff/oi = 0, (38a) 

foi qii - fqw- cigio qu + C 2 g qoi = 0, (38b) 

foiPoi- fp- 02910Pio + cigp = 0, (38c) 

foi- f - Cl (510 - 9 ) - f 910 Pio qio + 9 foi Poi qoi = 0. (38d) 

We can solve equations (l38l) for poi, qoi, foi and 9 (or for pio, qio, f and 510 )■ If ci = C 2 = 1, we 


derive two sets of solutions, as in the case of the nonlinear Schrodinger. Specifically, the first branch is 
the singular solution already given in (I17p . while the second branch involves rational expressions of the 
remaining variables. When either ci or C 2 is equal to 0, then system (1381) admits a unique non-trivial 
solution. This solution is given by 


POI 

qoi 


(/V^io + C2fp{9ioPioqio - 1) - clgiopio + ciC29iop) , 
^^^^2 (filii - qio + 9ioPioqioqii) + ci5iogii(5ioPiogii - 1)) 


/oi 
9 = 


= f 


B 
A' 
A 

910 


(39a) 

(39b) 


where A := fpqu + 02 ( 510^10911 - 1) and B := fpqio + cigiopqii - 02 - 

In this discrete context, the first integrals of the Backlund transformation given in (1311) become first 
integrals for system (l3^ . i.e. 


(T - 1) {fp qio - 02 f) =0, (5 - 1) {g^p qoi - 5 ) = 0. (40) 

Moreover, a generalised symmetry of the latter system follows from (|3Up and it is generated by the 
differential-difference equations 


dxP = p {pio qio -pq) - 

dxq = q{pq-p-io 9 - 10 ) 
dxf = f {pq-Pio 910 ), 


02 P 10 -oip , . poi - P 

- 7 - =p [Poiqoi - pq) -, 

/ 9 

ciq - C29-10 f \ ^0,-1 

= 9 (P 9 - 7 'o,-i 9 o,-i) - 


f -10 

dx9 = 9{pq-Poi 9oi)- 


9o,-i 


4.1.1 First integrals and a seven point scalar difference equation 

Let us consider now system (I38p with ci = C 2 = 1 and try to implement the first integrals (1401) so that to 
reduce the number of functions involved in this system by setting 

fpqio - f = a{n), 9^Pqoi - 9 = I3{m). (41) 

One option is to use the above relations to replace /, 9 in terms of p and q. In this case, we must solve 
equations (BH), which are quadratic in / and 9 , and hence introduce square roots, and finally derive a 
system of non-polynomial equations (correspondences) for p and q. 
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Another option is, instead of eliminating / and g, to use relations (1411) to replace the shifts of q. In 
this case, equations (f^T]) imply 


qio = 


a jn) + / 

Pp 


■) qoi — 


I3{m) +g 


g'^p 


Moreover, equation (I38ap suggests to introduce a potential u through the relations 


/ = 




g = 


Uoi 


(42) 


(43) 


Additionally, we introduce v by u := pju for convenience. 

Applying all the above substitutions to system ([38|) . we derive a system for v and u, namely 


uu + a{n)uoi uii + /3(m)uio 


= 0 , uii (uio - uoi) + V (uio - uoi) = 0, 


'^01 'WlO 

while a symmetry for this system is generated by 

—V U 'll 

dxU = -(u + a{n - l)u_io), dxV = -(u + a(n)uio) = -(u + /3{m)voi) . 

U_io UlO Uqi 


(44) 


(45) 


From equations ()44l) we can derive a higher order scalar equation either for u or for v, namely 
(5- 1) log A(n,u,u_io) - (T- l)logy(m,u,uo-i) + {ST~^ — l) logZ(n,m,u,u-n) = 0, (46) 


and 


(5- l)logX(n,u_io,u) - (T- 1) logy(m,uo-i,u) + {ST ^ - l)logZ{n,m,v-ii,v) =0, (47) 


where 


z — u 


X 'y 

X{n,u,x) = l + a{n-l)-, Y{m,u,y) = 1 +f5{m-l)-, Z{n,m,u,z) = — ^ ■ 

u u a{n — l)z — p(m)u 

A symmetry for equation (1461) follows from (|45D and it is generated by 

dxU = uX{n, u, u_io) Y{m + 1, uqi, u) Z{n, m, u, U-u) 


(48) 


while 

dxV = VX{n, V-io,v) Y{m + 1, v, uqi) Z{n, m, v-ii,v) 
generates a symmetry for equation (|47p . 

Equations (1461) . (|47l) are similar and have the same properties. They are defined on a stencil of seven 
points and can be solved uniquely with respect to any Uij and Vij except u and v, respectively. Because 
of this feature, if initial data are given along a double staircase, then these equations uniquely determine 
the evolution above and below this initial configuration as it is shown in Figure El 


Remark. When a{n), j3{m) are constants, i.e. a{n) = a, /3(m) = /3, equations (HSP and (ITTP are related 
to the discrete Toda equation 


(<S - 1) log(e"'-"'-i° +1) + (T - 1) log(e^-“-"' 


l) + (5r-l)log 


e"’-"’-i-i + 1 


= 0 , 


a 



(49) 
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The two quadrilaterals with one 
common vertex where equation 
is defined 


The initial value problem 


Figure 3: The stencil of seven points and the initial value problem 


i.e. equation (H) in [3]. This relation is made evident if we first reverse the m direction, i.e. change 
indices (i,j) to (i,—j) and operator T to its inverse in both equations (|16|) . (H71) . and then make 
the point transformation 


u = and v = (-l)”^a-^/3-”^e"', 


to each equation, respectively. In this context, system (IHl) defines the self-duality transformation for the 
Toda equation (j49p [3]. In particular, if we make the above change of variables to system ()44p . then it 
will become 




qUIIO—100,-1 


+ 7 


gUlio-UIo,_l _j_ 






- 1 ) 


^WQ^-l—W _ _ 


1 0^10—^ 
ry ^ 0,-1 _|_ \ 




1 + 1) , (50) 


where w and w are two different solutions of equation (I49p . 


4.1.2 Lax pair with matrix fl32p and a six point difference equation 

Let us consider now the Lax pair 


4'io = M(p)4', Toi = 

where matrix M{p) is given in (I32p and M{p,qQi,g]l,l) in (I29|) . It follows from the compatibility 
condition of the above pair that 


m =P, 9 = 


Poi-P P 

-T) 510 — - 9, 

p[poiP-ii - ppw) Pm 


and finally we arrive at the six point difference equation 

Pii - Pio Poi - P 


PioiPiiPoi - P 10 P 20 ) PoiiPoiP-ii - PPio) ’ 

We also find a first integral and a symmetry of this equation, which are given by 

(P -Poi)(Pio -P- 11 ) 


(5-1)- 




= 0 and d^p = p^ {pio - p_io) 


( 51 ) 


2 







































The stencil where the equation is defined 


The initial value problem 


Figure 4: The stencil of six points and the initial value problem for equation \51\) 


respectively. 

It is worth noting that equation (1^ can be uniquely solved with respect to any value of p except 
pio and poi- If initial data are given along a double staircase as it is shown in Figure 01 which must 
be consistent with the first integral, then the evolution of these data is uniquely determined above and 
below the double staircase by equation (ICT) . 

Remark. If we set the value of the above first integral to a{m) and, subsequently, make the change of 
independent variables {n,m) i—>■ {k,l) := (n + m,m), then we will arrive at the following quadrilateral 
equation for p{k, 1) = p{n, m). 

{P-Pii) (Pio-Poi) = a(0 {PPW - PoiPiif ■ 


5 1^2 reduction group: Generic orbit 

A Z 2 invariant Lax operator with simple poles in the generic orbit can be taken in the form 


C — Dx + 


I 


-S - 


1 


A - 1 A + 1 
The corresponding NLS type equation is 


o'aS'fTs, 


S := 


1 


p-q 


Pt - Pxx 


2 pl ^ Spqpx - ^p'^qx 


p-q 


{P - Q? 


Q.t — Qxx 


p + q —2pq 
2 -p-q 


2ql ^ 8qpqx - ‘iq^Px 


(52) 


p-q 


{P - <lf 


which is actually equation (m) in 

The Darboux matrix for the above Lax operator is derived in the same way as in the previous section 
and three distinct cases occur. 


1. The first Darboux matrix is 
0 oi 


M = A 


+ 


hiu 0 
0 b 2 V 


02 0 

and the Backlund transformation is given by 

ai + bipu ai + biqu 

Pw = T-;-) Qw = 


uv = 1, ai,bi€^C and |aia 2 |^ + 1 ^ 1 ^ 21 ^ 7 ^ 0, (53) 


b 2 V + a 2 P 


b2V + a2q 


drU = 4 u 


Pio 


P 


pio -qio p-q 


(54) 
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This transformation contains as particular subcases two Darboux transformations related to point 
symmetries, namely scalings (ai = 02 = 0 , n = n = 1 ) and inversions (61 = 62 = 0 , ai = 02 = 1 ). 

2. The second Darboux matrix is 


Mip,qio) = 


1 


A-1 

and the Backhmd transformation is given by 

-1 


Qio 1 

1 —p 


qio -1 

A +1 V -1 -p 


qio = —, dxP = 4p 

p \p-q Pio- qio 


PlO 


3 . The last Darboux matrix is given by 

/ f qio -pqio 


M{p,qioJ-,ci,C2) = 


A-1 


1 


-P 


qio pqio 


A + 1 V -1 -P 


+ 


Cl 0 
0 C2 


(55) 


(56) 


(57) 


where ci, C 2 are constants such that |cip + |c 2 p 7 ^ 0 and, without loss of generality, we can set 
these constants equal to 0 or 1. The derivatives of p, qiQ and / are given by the following relations 


Function 


dxP 


dxqio 


dxf 


4p 

-4 


Pio 


qio 


qio - pio 

PlO 


qio - PlO 


P \ ^ C2P10 - cip 

q-pj f qio-PlO 

_ P \ 2 C2gio -ciq 

q-p) f q-p 


2 ci 2 c2 


qio - PlO q-p 


4'(ci,C 2) = (2/p + C 2 ) (2/gio - Cl) 


(58) 


(59) 


defines a first integral for equations (1^51) . i.e. Dx^{ci,C 2 ) = 0 on solutions of the latter system. 


5.1 Derivation of discrete systems 

The first discrete Lax pair to consider is 

Tio = M{p, qio, f] Cl, 02)41 , Toi = M{p,qoi,g-,l,l)4>, (60) 

where matrix M is given in (1571) and for constants ci, C 2 one may consider three distinct cases : (i) 


Cl = C 2 = 1, (ii) Cl = 1, C 2 = 0 and (hi) ci = 0, C 2 = 1. 

In this generic setting, the compatibility condition of system (1601) results to 

hi- f - cigio + C 2 g = 0, (Ola) 

/ 01 P 01911 - f pqio- 02910 PlO qii +cigpqoi = 0, (61b) 

foiqii - fqio - Cigioqii + cigqoi - 2qii{ghiqoi - Piofqio) = 0, (61c) 

foiPoi - fp- 02910 Pio + C 2 gp + 2p{g hiPoi - 9 io f Pio) = 0. (61d) 


This system can be solved for either (poi; 9oi)/oi; <?) or {pio,qio, f, 9io)- When ci = C 2 = 1, then it leads 
to a solution with two branches: one branch is the trivial solution CB, while the non-trivial branch 
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involves rational expressions of the remaining variables. In the other two cases (ci = 1, C 2 = 0 or ci = 0, 
C 2 = 1), system ()61l) admits a unique non-trivial solution. In any case, the non-trivial branch can be 
easily found, but is omitted here because of its length, and we consider it as a difference system. For this 
system, it can be verihed that it admits two first integrals. 


(r-I)(2/p + C2)(2/^io-ci) =0, {S-l){2gp+l){2gqoi-l)=Q, (62) 

and a symmetry generated by 


= 2p 

dxQ = - 2 q 


Pio 


P 


qio -pio q-pj f qio -pio 


1 C2P10 - Cip 
+ --= 


Poi 


P 


qoi -poi q-p 


+ 


1 Poi -p 

9 qoi - Poi 


p 


P-IO 


q-p q-io-p-ioj f-io q-io-p- 


+ 


1 C2q - ciq- 


10 


-10 


(63) 


= -2q 


P 


Po-i 


dxf = 


Cl 


q-p qo-i-po -1 

C2 


+ 


1 q- qo,-i 


qio -Pio q-p 


■) dxg — 


9 o,-i qo-i -P-io 

1 1 


9oi -Poi q-p 


We can use the two hrst integrals (j62p to reduce the number of dependent variables involved in system 
(j6ip . In particular, we have two different options. The first option is to use the hrst integrals to remove 
function q from the system and a conservation law to replace / and g with a potential u, as we did in 
the previous section. The second option is to consider particular values for these integrals so that to 
eliminate / and g. These considerations are presented in the following two subsections. 


5.1.1 First integrals and a seven point scalar equation 

Let us consider the case ci = C 2 = 1 for system ([HT]) and its integrals (f62l) . Choosing the values of the 
latter, 

(2/p + 1 )( 2 /( 7 io- 1 ) = a(n)-l, {2gp + 1) {2g qoi - 1) = (3{m) - 1, (64) 

we can express gio and goi in terms of p, f and g as 


qio 


1 2 fp + a{n) 

V 2/p + l 


1 2 gp + f3{m) 
2 g 2gp + 1 


(65) 


Moreover, the hrst equation of (I6ip for ci = C 2 = 1 has the form of a conservation law, suggesting the 
introduction of a potential u via the relations 


f = uio-u, g = uoi -u. 


( 66 ) 


We use now relations ()65p . ()66l) to eliminate q, f and g from system (|6ip and derive the following system 
for p and u. 


2pio 

2poi 


a{n) — /3{m) j3{m) 2(a(n) — l)p 

uio - uoi till - uio 1 + 2p{uio - u) ’ 

a{n) — /3{m) a{n) 2{j3{m) — l)p 

uio - uoi till -uoi 1 + 2 p{uoi - u) ' 


(67a) 

(67b) 


A symmetry of this system easily follows from (l63P by using substitutions (l65P , (f66P but it is omitted here 
because of its length. Equations (l67P can be solved uniquely either for the pair {pio, uio) or for (poi; "^^oi)) 
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but here we present it in this form because it is more elegant and convenient. Moreover it makes apparent 
the invariance of the system under the involution {pij,Uij, a{n), f3{m)) o {pji, uji, f3{m),a{n)). Regarding 
the Cauchy problem, initial values along a staircase are compatible with the evolution defined by the 
above system. 

Equations (l67|l can be decoupled to a scalar equation for u. Indeed, the compatibility condition 
T (Pio) = S (poi) implies that u must obey the equation 


(5 - 1) - (T - 1) + (ST-' - 1) - 1) = 0, 


u - U-io 


U — Uo,-l 


U — U-ll 


( 68 ) 


which, up to point transformations, is the non-autonomous version of the Toda-type equation (A) in [3], 
cf. also [1]. A symmetry of this equations follows from the symmetry of system (1671) and is generated by 


dxU = 


{uiQ - u){uq-i - u){ui-i - u) _ (u_io - u)(uoi - - u) 


Ff 


0,-1 


F. 


-10 


(69) 


where 

Foo := a{n){u - moi)(iiio - mi) - I3{m){u - uio){uoi - un). (70) 

Equation (l68l) is defined on a stencil of seven points and can be solved uniquely with respect to any 
Uij except u. Because of this property, if initial data are given along a double staircase, then equation 
(j68l) uniquely determine the evolution above and below this initial configuration as it is shown in EigureEl 


Remark. Equation (|68p is the Euler-Lagrange equation for the Lagrangian 

C = Q:(n - 1) log(n - n_io) - /3(m - 1) log(u - uq- i) - (a(n - 1) - /3{m - 1)) log(u_io - uo-i), 

which is also considered as a Lagrangian for the discrete Schwarzian KdV or QIq |18] . the form of which 
is Foo = 0 pale], where F is given in (I7np . 


5.1.2 First integrals and a five point scalar equation 

Now we consider the case ci = C 2 = 1 and two particular values for the first integrals given in (|62p . More 
precisely, let us consider that 


(2/p + l)(2/gio - 1) = 0, {2gp + l){2gqoi - 1) = -1, (71) 


from which we can express / and g in terms of p and q rationally. While the second equation determines 


g unique! 


3, the first equation admits two different solutions and we choos^ 


-1 _ 1 M 

2p ’ ^ 2 Vgoi pj ' 


(72) 


Then, for ci = C 2 = 1 and in view of substitutions ([72]) . system (I6TI) and its symmetry (f63P reduce to 


Pi 0 -P = 910 - 901,-=- (73a) 

Pio Poi 911 901 

^Tlie solution ^ = 0 is not considered since along it system m and its symmetry (ESP degenerate. 

^The second choice for / leads to a system related to (|73P by a point transformation. 
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and 


dxP = 


910 - Pio 


P 


dxQ 


pq 

p-q 


P-ioQ-io 
P-io - 9-10 ’ 


(73b) 


respectively. 

It can be readily verified that the above discrete system for p and q can be written in a conserved 
form as 

(S-l)(,-p) = (T-l),, (S-l)(i-J-) =(T- 1)1 

\P 901 / P 

We can use either of these conserved forms to introduce a potential and then derive an equation only 
for the potential. In either of the cases, we end up actually with the same scalar equation. Here, we 
introduce potential w employing the first conservation law and, in particular, we set 


p = wo-i-w_i^o, q = wo_i-w_i_i. 


The substitution of the above expressions into equations (1731) results to a scalar equation for potential 

tcjl 

ii^-, (74) 

w — icio w — in_io w — wi-i w — in_ii 

and a symmetry of this equation is generated by 


dxW = 


{w - W-io){w - W-u) 

W-io - W-ii 


A staircase initial value problem for equation (1741) is similar to the one we considered for the Toda equa¬ 
tion in the previous section. That is, initial data can be given at points Wi-i and Wi-i-i from which a 
solution can be uniquely determined on the whole lattice. 


Remark. By the change of independent variables (n,m) {k,l) '■= [n + m,m), equation ()74l) can be 

written as “the missing identity of Frobenius” for the function w{k,l) = w{n,m), 


1 1 

:-1—:-:- 

W — WlQ W — W-IQ 


1 1 
:- h :-, 

W — Wol W — U)o,-l 


which appears in the theory of Fade approximants m, as well as in relation with the discrete KdV 
equations HI, H3 [Ml E] and the e-algorithm [27| . 


5.1.3 A Lax pair with matrix (I55|l and a six point difference eqnation 

Now we consider the discrete Lax pair 

Tio = M(p,gio)T, Toi = M{p,qoi,g;l,l)^, 

where M{p,qio) is given in (l55|) and M{p,qQi, g; 1,1) in (l571) . The compatibility condition of this system 
implies 

^ ^ ^ P-ii(Poi -p){'^+PPio) ^ {Poi -P)(l +F01P-11) 

p’ ^ 2 p(ppio-poiF-ii) ’ 2 (ppio-poiF-11) 

®If we used the second conservation law to introduce the potential, then the resulting equation would be related to (|74l) 
by interchanging n and m, i.e. changing indices (ij) to {ji). 
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which subsequently leads to the scalar difference equation 

PioPoi I {Pn +P-ii){pPioP 20 +Poi) - {p+P 2 o){P-iiPoiPii +Pio)} + (1 -PioPoi){pPioP 20 -P-iiPoiPii^ = 0. 

(75) 

This equation admits the first integral 


, ppiopoip-ii / I 1 

9 := - -^ pH-poi- 

(ppio -poiP-ll)^ V Pio P-11 


1 1 

- +P10-P-11 

P Poi 


(76) 


Moreover, a generalised symmetry of (j75p is generated by 

dxP = P ' 


1 


1+ppio l+pp_io. 

Finally, it can be easily shown that a non-autonomous symmetry of equation (I75p is generated by 


drP = 


n 


n — 1 


1 + PPio 1 + PP- 


10 




Equation ([Tsp is defined on a stencil of six points, cf. Figure [H and can be uniquely solved with 
respect to any value of p except pio and poi- Initial data for equation (ITHp can be given along a double 
staircase as it is shown in Figure 01 

Remark. If we set the value of the first integral (I76p to a(m) and, subsequently, make the change of 
independent variables {n,m) i—>■ {k,l) := {n + m,m), then equation = a(m) will become a quadrilateral 
equation (correspondence) for p{k,l) = p(n,m), namely 


ppioPoiPiiH{p)H [p = a{l){ppio -poipnf 


H{p) :=p + ^ ^ - pn- (77) 

Pio Poi 


Obviously if we set a{l) = 0, the above equation reduces to Hirota’s discrete KdV equation [T2] either 
in the form H{p) = 0 or H{p~^) = 0. Hence, we consider equation (1771) as a quadratic Hirota KdV 
equation. This relation allowed us to derive the non-autonomous symmetry of equation (I75h from the 
corresponding symmetries of Hirota’s KdV equation m- 


6 Dihedral reduction group: Degenerate orbit 

We now consider Lax operators which are invariant with respect to the following transformations 

si(A) : T(A) ^ fj 3 /:(-A)fj 3 , S 2 (A) : T(A)uiT (A“^) cji, ui = ^ ^ (78) 

Here, the reduction group is generated by the above set of involutions and it is isomorphic to Z 2 XZ 2 — B 2 . 
The invariant Lax operator corresponding to the degenerate orbit can be taken in the form 


C - Dx+ + A ( 


1/0 2q\ 1 

+ A I 2p 0 J ■ 


(79) 


This operator corresponds to the following deformation of the derivative NLS equation 

Pt=Pxx + 8{p^q)x-4.qx, qt =-Qxx + 8{pq^)x - 4.px. (80) 
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It is simple to check that matrix <73 is a Darboux matrix for operator (1791) and corresponds to the 
discrete symmetry (p, q) !->■ (— p, —q) of system (IHOl) . A A-dependent Darbonx matrix for operator ([7^ is 


= “ (( 0 a “0 ^ K <110 0 ) ^ 

and the corresponding Backlund transformation is given by 


1 0 \ 1 / 0 qio\\ 

0 I X\p 0 JJ 


(81) 


dxP = ‘^(^{Pioqio-pq)p+{p-Pio)f+ q-Qio), (82a) 

<9x^10 = ‘^(^{pioqio-pq)qio + p-pio + {q-qio)fy (82b) 

dxf = ‘^(^ipwqw-pq)f+ ip-pio)p + iq-qw)qioy (82c) 

dxU = -2{pioqio - pq)u. (82d) 

It is straightforward to show that these differential equations admit two first integrals dx^^^^ = 0, i = 1,2, 
where 

$(i) =u^{f-p , $( 2 ) =u‘^ ^ (83) 

which imply that matrix M has constant determinant since 


det M 



(^(1) _g <|)(2)^ 


6.1 Derivation of discrete systems 

We introduce the discrete Lax pair 


^-10 = M(p,gio,/;M)^', -^01 = M{p,qoi,g-,v)^, (84) 

where matrix M is given in (1811) . The compatibility condition of this Lax pair leads to a set of equations 
for p, q, f and g, 


f 01 - f - 910 + 9+ Poi qoi - Pio qio = 0, (85a) 

(/oi - 9io)p + 9P01 - fPio + qoi - qio = 0, (85b) 

if - 9)qii + 910 qio - foi qoi - Poi + Pio = O , (85c) 

foi9 - f 910 + P {poi - Pio) + qii {qoi - qio) = O, (85d) 

and an equation solely for u and v, 

uoiv — viou = 0 . ( 86 ) 


Functions u, v are apparently redundant since they are completely separated from the remaining ones 
and are involved only in equation () 86 p . Taking the value of the first integral in ()83p to be 1, then we 
can set 

2 1 2 1 
u = - -, V = -. 

/ - pqio 9 - pqoi 

In view of this substitution, equation ([ 86 P becomes 

(T-l)ln(/-pgio) = (5-l)ln(5f-pgoi), 
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which can be easily verified to be a conservation law for equations (1851) . 

Equations ([851) can be easily solved with respect to (poi> 9oi>/oi> ff) or {pio, qio, f, gio) leading to a 
solution with two branches: the trivial branch (I17p and the non-trivial one which we consider as a system 
of difference equations. For the latter system it can be easily verified that it admits two first integrals 


ir-1] 


f - pqw 


P - (p2 + g2^) + 1 


= 0 , 


(5-1)- 

9^ 


9 - pqoi 
{p^ + qoi) +1 


a conservation law 

{T - l){f +pq) = {S -l){ 9 +pq), 

and a symmetry given by 


dxP = {pioqio-pq)p+{p-Pio)f+ q-qio, 
dxq = {pq-p-ioq-io)q + P-io-p+{q-io-q)f-io, 
dxf = {pioqio-pq)f+ {p-pio)p + {q-qio)qio, 
dx9 = {poiqoi-pq)9 + {p-poi)p+{q-qoi)qoi- 


(87) 


Now, we will consider two particular values for the first integrals (1871) which allow us to reduce the 
number of functions involved in system (j85j) by expressing /, 9 polynomially in terms of p and q. 


6.1.1 First reduction and a Toda type eqnation 

Let us first consider for the first integrals the values 

/ - pqio _ ^ 9 - pqoi _ 1 

P - (p2 + g 2^) + 1 ’ ^2 _ 0^2 g 2 ’ 

which imply that 

f - pqio = 0, (5-p + goi - l)( 5 +p- 9 oi - 1) = 0. (88) 

From these algebraic equations, we choose the solution 


/ = pqw, 9 = p- qoi + i- 


(89) 


If we substitute these expressions into system ()85p . its conservation laws and symmetry and then make 
the point transformation {p,q) = {p — l,q — 1), we will come up with the system 


Poi 


Pwqw 

qii 


qoi 


(p-2)(gio -2)qn 

Pwqw — 2gii 


(90) 


along with its conservation laws 


(T - l){p - l)(gio + q- 2 ) = {S - l){pp - 1) - goi), (T - 1) Inpqio = {S - 1 ) Inp 


and its symmetry 


dxP = p{p - 2 ){qio - q), d^q = q{q - 2 ){p - p-w). 


Remark. Using the second conservation law above to introduce a potential w by 

p = exp {w - rco- 1 ), q = exp {wq-i - W-io ), 
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we derive the scalar equation 


^Woi-W _ gUI-UIO,_l _ 


1 

2 




— e 




and its symmetry 


dxW 




2 


(91) 


6.1.2 Second reduction and a seven point scalar equation 

Another choice for the values of the first integrals (j87p is 

f - pqio _ -1 g - pqoi _ i 

f2 _ (^p2 2 ’ 5,2 _ (p 2 + q2^-^ +1 2 ’ 

or, equivalently, 

if+p + qio + l){f-p-qio + 1) = 0, {g-p + qoi-l)ig+p-qoi-'i-) = 0. (92) 

The above equation has four solutions of solutions and we choose 


/= P + 910-1, 5 =P- 901 + 1- (93) 

As before, the substitution of (|93]l into equations (1851) and the point transformation {p, q) i-)- (p + 1 , g + 1 ) 
result to 


Poi = Pio - 911 + 2 + 


Pio(9io - 2) 


P 


901 = P 


Pio9io — 2911 


Pio(p + 9io - 2) -P9ii’ 


(94) 


Similarly, we find two conservation laws 


(T - l)pq = {S - l){pq - 2goi), (T - 1) ln(p - 2)(gio - 2) = (5 - 1) lnp(goi - 2), (95) 

and a symmetrjQ 


dxP ={p- 2)(pio(9io - 2) -p(g - 2 )), dxq = {q - 2){q{p - 2 ) - g_io(p-io - 2)). (96) 


Remark. Using the first conservation law in (I95p . we can introduce a potential y by 

p = 2 ——, g = 2 ( 2 / 0-1 - ?/_io), (97) 

yo-i - y-10 

and derive the Toda-type equatiorU 

(<S - 1) log ( 2 / - 2 /- 10 ) - (T - 1) log ( 2 / 0-1 -y) + {ST~^ - 1 ) log Tl- - -^ = 0 (98) 

V y-y-iij 

^These differential-difference equations are related to the relativistic Volterra lattice [30] 

dyP = P (Qio — Q + hPioQio — hPQ ), dyQ — Q{P — P_io -|- hPQ — hP-ioQ- 10 ) 

by the point transformation {x,pij,qij) e-l iy,Pij,Qij) ~ ^^^ 2 ^’%^)• 

®If we reverse the m direction, then this equation will become equation (D) in [5], cf. also [?]. 
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A conserved form of equation (|98l) is given by 

(5 _ 1) In (^0-1 ~ y ) (^1-1 - 2/ - 1) ^ (-7- _ ;l) In - y + 1 ) 


U-w — Vo-i 

while the differential-difference equation 


y-10 - 2/0,-1 


dxv = (yio - y)(yo,-i - y) ( — -1) = -(2/01 - y){y-io - y) (— -h 1) 

\yi-i-y J \y-ii-y J 

dehnes a symmetry of this equation. 

Remark. If we use the second conservation law in ()95p to introduce a potential (j) by the relations 

2 


P = 


X _|_ e<2>o,-i-</'-io ’ 
then system (fM)l will reduce to equatior]^ 


q 


= 2 (1 + ^ 


-10—(p-1-1 


(99) 


(5 - 1) e'^ 


+ (r-l)e' 




( 57-^1 - 1 ) 


1 


A symmetry of this equation is generated by 


dU = + 


1 


1 -h 

= X -(- _ g<2>io-<2> _ 


= 0 . 


( 100 ) 


1 


Remark. Combining transformations (I97p and (I99p . we derive the duality transformation [3] 

1 1 


2/10 “ 2/ = + - 7 . -J 

which connect solutions of equations (fMP . (llOOp . 


, 2/01 - 2/ = -e" 


1 + e'P-i.i 


( 101 ) 


7 Concluding remarks 

In this paper we discussed the Darboux-Lax scheme for Lax operators related to nonlinear NLS type 
equations. We derived integrable systems of differential-difference and partial difference equations, and 
discussed several reductions to scalar Toda-type equations. The results of this paper have already been 
employed in other works. For some of the systems presented here the symmetry structure and recursion 
operators were studied in m , whereas the connection of our systems with Yang-Baxter maps was explored 
in m- Lax-Darboux schemes corresponding to other Lie algebras are discussed in [11122]. 
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®If we reverse the m direction, the resulting equation will be equation (E) in [3] . 


23 















References 

[1] M. Ablowitz and H. Segur (1981), Solitons and the Inverse Scattering Transform, SIAM, Philadelphia 

[2] V. E. Adler (1998) Backlund Transformation for the Krichever-Novikov equation Int. Math. Res. 
Notices 1998 1-4 

[3] V. E. Adler (2000) On the Structure of the Backlund Transformations for the Relativistic Lattices 
,J. N. Math. Phys. 7 34-56 

[4] V. E. Adler (2001) Discrete equations on planar graphs J. Phys. A : Math. Gen. 34 10453 

[5] V. E. Adler (2010) Classification of discrete integrable equations D.Sc. Dissertation, Russian 
Academy of Sciences, Institute of Theoretical Physics L. D. Landau (in Russian) 

[6] V. E. Adler, A.I. Bobenko and Yu.B. Suris (2003) Classification of integrable equations on quad- 
graphs. The consistency approach Comm. Math. Phys. 233 513-543 

[7] V.E. Adler and R.I. Yamilov (1994) Explicit auto-transformations of integrable chains J. Phys. A: 
Math. Gen. 27 477-492 

[8] G. Berkeley, A. Mikhailov, P. Xenitidis (2014) in preparation 

[9] R.T. Bury (2010) Automorphic Lie Algebras, Corresponding Integrable Systems and their Soliton 
Solutions, PhD thesis. University of Leeds 

[10] W. B. Gragg (1972) The Pade Table and Its Relation to Gertain Algorithms of Numerical Analysis 
SIAM Review 14 1-62 

[11] J. Hietarinta and G. Viallet (2011) Integrable lattice equations with vertex and bond variables J. 
Phys. A Math. Theor. 44 385201. 

[12] R. Hirota (1977) Nonlinear Partial Difference Equations. 1. A Difference Analogue of the Korteweg-de 
Vries Equation J Phys Soc Japan 43 1424-1433 

[13] D.J. Kaup and A.G. Newell (1978) An exact solution for a derivative nonlinear Schrodinger equation. 
J. Math. Phys. 19 798-801 

[14] F. Khanizadeh, A.V. Mikhailov and J.P. Wang (2013) Darboux transformations and recursion op¬ 
erators for differential-difference equations. Theoret. and Math. Phys. 177 1606-1654 

[15] S. Konstantinou-Rizos and A. V. Mikhailov (2013) Darboux transformations, finite reduction groups 
and related Yang-Baxter maps J. Phys. A: Math. Theor. 46 425201 

[16] D. Levi (1981) Nonlinear differential difference equations as Backlund transformations J. Phys. A: 
Math. Gen. 14 1083-1098 

[17] D. Levi and R. Benguria (1980) Backlund transformations and nonlinear differential difference equa¬ 
tions Proc. Natl. Acad. Sci. USA 77 5025-5027 

[18] S. B. Lobb and F. W. Nijhoff (2009) Lagrangian multiforms and multidimensional consistency J. 
Phys. A: Math. Theor. 42 454013 (18pp) 


24 



[19] S. Lombardo and A.V. Mikhailov (2005) Reduction groups and automorphic Lie algebras Commun. 
Math. Phys. 258 179-202 

[20] S. Lombardo and J. Sanders (2010) On the Classification of Automorphic Lie Algebras Commun. 
Math. Phys. 299 793-824 

[21] A. V. Mikhailov (1981) The reduction problem and the inverse scattering method Physica D 3 
73-117 

[22] A.V. Mikhailov, G. Papamikos, and J.P. Wang (2014) Darboux transformation with Dihedral reduc¬ 
tion group J. Math. Phys. 55 113507 

[23] A.V. Mikhailov, A. B. Shabat, and R. 1. Yamilov (1988) Extension of the Module of Invertible 
Transformations. Classification of Integrable Systems Commun. Math. Phys. 115 1-19 

[24] F.W. Nijhoff and H.W. Capel (1995) Backlund Transformations and the Linearization of the 
Anisotropic Heisenberg Spin Chain Phys. Lett. A 91 431-434 

[25] F. W. Nijhoff and A. J. Walker (2001) The discrete and continuous Painleve hierarchy and the 
Gamier system Glasgow Math. J. 43A 109-123 

[26] S. P. Novikov, S. V. Manakov, L. P. Pitaevskii and V. E. Zakharov (1984), The theory of solitons. 
The inverse seattering methods, Plenum Publishing Corporation. 

[27] V. Papageorgiou, B. Grammaticos and A. Ramani (1993) Integrable lattices and convergence accel¬ 
eration algorithms Phys. Lett. A 179 111-115 

[28] G.R.W. Quispel, F.W. Nijhoff, H.W. Capel and J. van der Linden (1984) Linear integral equations 
and nonlinear difference-difference equations Physica A 125 344-380 

[29] Yu. B. Suris (1995) Bi-Hamiltonian structure of the qd algorithm and new discretizations of the 
Toda lattice Phys. Lett. A 206 153-161 

[30] Yu. B. Suris and Or. Ragnisco (1999) What is the Relativistic Volterra Lattice? Comm. Math. Phys. 
200 445-485 

[31] P. Xenitidis (2012) Linearization and Potentiation of difference equations, Talk in SIDEIO (Symme¬ 
tries & Integrability of Differenee Equations), Xikou, Ningbo, China, 11-15/06/2012 

[32] V.E. Zaharov and A.B. Shabat (1972) Exact theory of two-dimensional self-focusing and one¬ 
dimensional self-modulation of waves in nonlinear media JETP 34 62-69 


25 



